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I. INTRODUCTION 

Quantum field theories (QFT) formulated in non-commutative spaces have interesting phenomenological impli¬ 
cations which have been the subject of numerous investigations in recent years, in connection with many physics 
phenomena related to the possibility of violation of Lorentz and CPT symmetry [1-4], analyzes of matter-antimatter 
asymmetry [5] , implications in cosmological models [6-9] , in addition to observations involving ultra-energetic cosmic 
rays, neutrino physics, and other sources of high energy existing in the Universe [10, 11]. 

Such theories, also known as non-commutative quantum field theories (NCQFT), were inspired by quantum me¬ 
chanics on a non-commutative spacetime [12, 13] whose proposal is based on the idea of quantization of spacetime 
through the replacement of its coordinates by operators which meet the commutation relation This 

causes the spatial coordinates to become discrete and there is no longer the notion of points in the space. This 
proposal follows the same lines of reasoning in the usual quantum theory characterized by quantization of the phase 
space which results in the existence of a minimal area scaled by Planck constant. 

The NCQFT was initially motivated by the interest of eliminating divergences that appear in quantum field theory, 
since by introducing non-commutativity is possible to exist a minimal distance in the spacetime to eliminate the 
infinities [14, 15]. However, this proposal was abandoned because of the success of renormalization process. 

In addition, the non-commutative Yang-Mills theory emerged later as a limit of low energy of string theory in the 
presence of a magnetic field background [16], causing a major boost in the resumption of the study of models involving 
non-commutative QFT. The reappearance of these theories is also motivated by some formulations of quantum gravity, 
which suggests that in regions around the Planck scale (Ip ^ 10“^^cto) the spacetime loses its continuous structure 
and the effects of quantum gravity become important [17-19]. 

An important aspect presented by non-commutative quantum field theory is the violation of the Lorentz invariance, 
which has been the subject of great attention in QFT [1, 2, 20]. This violation should occur only at very small distances, 
thus, within the limits that the non-commutative parameter tends to zero the results previously known are recovered. 
Theories with loss of Lorentz invariance present deformations in the dispersion relation. These deformations allow 
us to build phenomenological models that can be further tested in cosmology or particle physics. In this perspective, 
it is expected that signs of the non-commutativity appear in experiments involving cosmic microwave background, 
ultra-energetic cosmic rays, or other source of high energy. 

The introduction of the non-commutativity in QFT can be made at the fields by modifying their canonical com¬ 
mutation relations [1, 5, 21], which we will discuss later in the section 11. These fields, called non-commutative fields, 
arise as a generalization of quantum mechanics in Groenewold-Moyal plane [22], with the target space (or space of 
fields) being considered as a non-commutative plane and a base space considered a commutative spacetime. 

Several authors investigated the physical implications caused by these non-commutative fields. In Ref. [1] was 
shown how the non-commutativity in the field space produces violation of Lorentz invariance which is manifest in the 
dispersion relation. Subsequently, in Ref. [23] was shown how the blackbody radiation spectrum is modified by the 
deformation of the canonical commutation relations of massless scalar fields. An important aspect of this modified 
spectrum is that the energy density in the region of high frequency is higher than that of the usual radiation blackbody 
spectrum. Later, in Ref. [6] was presented implications of this theory in inflationary cosmological models in a universe 
filled with a non-commutative gas at high temperature regime. Also, in this perspective, recently in Ref. [24] has been 
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studied how the Casimir effect might be affected by the non-commutativity in the target space, where was shown the 
arising of a repulsive non-commutative Casimir force at the microscopic level. 

In this paper, we consider this approach based on non-commutative fields with commutation relations in equal 
times given by 

= ie“^6»(cr;f-y) (1) 

to investigate the Bose-Einstein condensation of a system described by non-commutative scalar fields. Here 0 ((t, x — y) 
is considered as a Gaussian distribution with a parameter a related to the standard deviation of the distribution. 

The Bose-Einstein condensation has some very peculiar characteristics, such as below certain temperature a finite 
fraction of the number of particles are condensed in the ground state. Theoretical and experimental studies discussing 
the properties of this thermodynamic phenomenon has gained great attention in different areas of physics [25-29]. 

Recently several studies have been developed in the analysis of relativistic effects in Bose-Einstein condensates [BO¬ 
SS] — see also [34] for related issues. However, it is also interesting to explore possible non-commutative consequences 
that become relevant for these systems at high energies. Because of this, we study the thermodynamic properties of 
bosonic fields described in the non-commutative quantum field theory assuming the inclusion of particle-antiparticle 
pair production. 

The paper is organized as follows. In Sec. H we make a revision on field theories in non-commutative spaces 
involving massive scalar fields, where the physical peculiarities of such theories based on the non-commutative algebra 
are analyzed. In Sec. HI we explore the effects due to the non-commutativity in the thermodynamic functions of 
a relativistic bosonic gas with and without anti-bosons. Expressions for these quantities were obtained, such as, 
particles density, internal energy, pressure and specific heat. In Sec. IV we present our results obtained through 
non-commutative theory which are compared to the results established by the usual theory. Finally, in Sec. V we 
present our final comments and conclusions. 

II. NON-COMMUTATIVE QUANTUM FIELD THEORY 

Field theories based on concept of non-commutative target space is gaining more attention and interest as matter of 
study, whose properties and implications caused are studied and applied to different physics problems [1, 6, 21, 23, 35]. 
The formalism that describes non-commutative fields was proposed from an analogy with the fundamental ideas of 
the non-commutative quantum mechanics (NCQM), which suggests the replacement of non-commutative relationships 
between the coordinates of a spacetime to non-commutative relationships between the fields. 

The aim of this section is to present a revision of the non-commutative quantum field theory. The concepts here 
formulated will be fundamental to study the thermodynamic properties of a bosonic gas with its dispersion relation 
modified by the introduction of new parameters related to the non-commutativity. This section is divided into two 
subsections. The first subsection focus on the definition of non-commutative spacetime and the following subsection 
discusses an extension of this theory to a massive scalar field where the target space is considered a non-commutative 
space living in a commutative spacetime in (3 -I- l)-dimensions. 
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A. Non-commutative spacetime 


The notion of a non-commutative spacetime had already been proposed since the beginning of the quantum mechan¬ 
ics [15]. Motivated by the extension of the usual commutation relations between position and momentum, the NCQM 
imposes non-commutation relations between the coordinates of the spacetime. The coordinates are substituted by 
Hermitian operators of a non-commutative algebra [40], satisfying deformed commutation relations 


[^^r] 

= ief^‘'9 = i9>^'', 

( 2 a) 

[Pu^pA 

= 0, 

( 2 b) 

[x^,pA 

= *50, 

( 2 c) 


where 9 is known as non-commutative parameter having dimension of area, and is a real anti-symmetric matrix. 

The non-commutativity of the spacetime causes effects in their geometrical structure, since the spacetime ceases 
to be continuous and becomes discrete, leaving aside the notion of point that is replaced by a minimal area, called 
Planck cell. This implies a set of uncertainty relations among its coordinates 

> ^[^‘'1. (3) 

Note that in the limit 0 —>■ 0 the usual commutative theory is obtained. 

In [23], see also [37-39], the non-commutative coordinates are related to the usual coordinates by the deformation 

x^ = x^- (4a) 

(4b) 


The deformation introduced above is applied to the coordinates, which are the degrees of freedom of the system. 
Similarly, in quantum field theory the degrees of freedom are the fields defined at every point in the spacetime, so we 
will introduce the non-commutativity into the fields which will engender the non-commutative target space. 


B. Free massive bosons description 

In this section, mainly following [23], we review the theory of commutative and non-commutative massive scalar 
fields and then we analyze the effects caused by the non-commutativity of the helds. 

1. Commutative scalar field 

We define here the commutative target space which will be important to compare with the effects that the non- 
commutative fields develop on the dispersion relation. For this, consider a theory with two scalar fields in a (3 -|- 
l)-dimensional base space and target space as a commutative plane then 

ip : M 3 X M — y 

{x,t) I- p{x,t), 
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where the field components are denoted by with i = 1,2. Let us assume that each spatial direction compactified 
in with radius R. This causes the field components to be periodic in the spatial coordinates {x,y,z), 

If" {x + R,y + R, z + R,t) = f" {x, y, z, t) = f" (x, t). (5) 

The compactification of space allows us write the components of the field f" {x, t) as a Fourier series 

( 6 ) 

n 

where n = ( 711 , 712 , 713 ), with rii € Z, and then the Fourier components are 

^n{t) = j^jd^xe-^-^--ip"{x,t). ( 7 ) 

The Lagrangian density describing massive scalars fields is given by 

£ = I ^ , 

i 

= ( 8 ) 

i 

whose Lagrangian is L = f C(Px 

^ = IE / ■ ( 9 ) 

i 

To obtain the above Lagrangian written in terms of Fourier modes we use the components of the fields f" (x, t) given 
in ( 6 ), then 

J (fx{dtip"f = 

n 

J d^x{Vf"f = £>-n, 

fi ^ 

J d^x{f"f = '^R^ip\if"_fi, 

n 

( 10 ) 


resulting in 




( 27 r|n| 


777 




( 11 ) 


The canonical momenta associated with the Fourier modes ip"z are 

, dL 

-^n = 

dft 

_ d .j .j 

2 ^ ^ m) 

j,rh 

2 / ^ r—m ' rm i n / 

= gR^v-n- 


(12) 
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Now using the Lagrangian (11) and applying a Legendre transform, we obtain the Hamiltonian H as a function of 
the modes tpb and canonical momenta ttI 


H = 

i,n 

f 

= E 


2gi?3 2 


27r|n| 

R 




= E 


2gi?3 


2 t i 

-TT^n^Pn^P-f 


(13) 


This Hamiltonian is equivalent to an infinite set of uncoupled harmonic oscillators with frequencies 


27r|n| 

R 


+ m?. 


(14) 


2. Non-commutative scalar field 


In this subsection, the commutative plane defined above is substituted by a non-commutative plane in . The way 
to introduce non-commutativity in the fields space is analogous to the deformation made on the spacetime coordinates 
presented in section (HA). The scalar field with non-commutative target space K^, denoted by t), can be written 
in terms of commutative fields as [23] 

= (/?“(f,t) - ^e“'’07rh(x,t), (15a) 

TTa{x,t) = 7Ta{x,t). (15b) 


Thus, the commutation relations in equal times for non-commutative scalar fields introduced above become (see 
appendix A) 


[ip‘^{x,t),(fi’’{y,t)] 

= ie°^^95{x — y), 

(16a) 

[Kaix,t),TTt,{y,t)] 

= 0, 

(16b) 

[(p°‘{x,t),nb{y,t)] 

= i5l5{x-y). 

(16c) 


In this scenario 0 is the non-commutative parameter with dimension of length. 

The function 95(x — y) that appears in Eq. (16a) has been regularized by a Gaussian-type distribution [23], written 
in terms of a new parameter a associated with the standard deviation of the distribution 
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0{(j) = 


{y/^a} 


-exp 


-E 


{Xj - yif 
2(t2 


(17) 


where 0{(j) = 9{a\x — y) = 95{x — y) in the limit cr —>■ 0, 9{a) has dimension of (length)”^, and we assume the 
simplification ai = a 2 = cfs = cr for the argument of the exponential. The behaviour of this function 9{a) can still be 
seen in the Figure [1] for different values of cr. 

We note that the regularization parameter a controls the width of the Gaussian describing the range of the non¬ 
commutativity, that is, in the limit as this parameter goes to zero we get the usual Dirac delta function and in the 
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0(a) 



FIG. 1: The function 6{a) for different values of cr, with 0 = 1.0. 


limit as the parameter grows we recover the commutativity of the helds. This regularization proved to be useful in 
study of the radiation spectrum of deformed black body [23] , where it was predicted a divergence in the energy density 
with respect to frequency increase when cr = 0 and that the deviation of the new radiation spectrum of the deformed 
blackbody has a stronger dependence on tr than 9. Besides, 9{a) can be written in terms of 

9{n) = 9e . (18) 


Analogous to Eq. (6) the non-commutative fields (,3“ (a;, t) can be written in terms of Fourier modes 




Then, using equations (15) we can get dressed transformations 

1 


‘Pfi = Pf, 
ttS = ttS 


2^3e“'>0(n)7r^^, 


(19) 


(20a) 

(20b) 


The commutation relations for the Fourier components become (appendix A) 


L'‘n5 '‘mj 


^3 ^n+m,0; 

0 , 

iS‘^^Sn,A- 


(21a) 

(21b) 

(21c) 


The non-commutativity of the fields introduced in (16) is equivalent to replace the Hamiltonian operator (13) for a 
new non-commutative Hamiltonian 


hnc = 


2gi?3 


9^^ 2 -i 
-^^nPnP-f 


( 22 ) 











Using the equations (20) we find 


hnc = 


2 ? ? 9 2 rt/ \ 4 k 

+ —UJaVn^P-n - 2^nQin)Uk^^'Kf, 


where 


Ui = 1 


gujft0{n) 


and Wf? = 


^27r|n| 


+ m^. 


(23) 


(24) 


Note from Eq. (23) that the dispersion relation for these boson fields was modified. If we take the limit 0 —0 we see 
that 9[n) —>• 0 and f2| = 1, falling into the usual dispersion relation, Eq. (13). 

It is remarkable that the Hamiltonian (23) has an equivalent term to a set of infinite harmonic oscillators and 
another proportional to the z-component of the angular momentum, i.e.. 


Hn = 


2 i , 

—^TT-TT* - H-- 

2,gR^ ^ 2 


T 2 \ '' ^ , A 

~ / . ^ikVn^n 

i,k 


(harmonic oscillators), 
(angular momentum). 


Thus, 


Hnc = E 


In addition, Hfi given by (25a) can be written in the usual form of the Hamiltonian of a harmonic oscillator 




and M = 


gR^ 


n 


2 ■ 


Hn = J2 

i 

with frequency u)ft and mass M given by 

UJfl — 

The hermiticity of the fields yi® implies that 

‘Pnit) = (^ / 

= gR^'Pn = t^-u- 


(25a) 

(25b) 

(26) 

(27) 

(28) 

(29a) 

(29b) 


So we can diagonalize Hfi by the introduction of creation and annihilation operators 



where 


As = Mujn = 


gR^ujn 


The commutation relations for 


and ab are (appendix A) 







0 , 

0 , 


(30a) 

(30b) 

(31) 

(32a) 

(32b) 

(32c) 
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The fields in terms of creation and annihilation operators, Eq. (30), become 



Substituting Eqs. (33) in Eq. (25a) we get Hfi in terms of creation and annihilation operators 



i 


(33a) 

(33b) 


(34) 


We can still write the operators in the normal order (represented by the symbol : :) by rearranging the creation and 
annihilation operators through the use of the algebra, that is 


E _ VlfiUjfi /lit, i t I \ 

^fi — / . 2 ■ 4 " J ■ 

fi i,n 



i,n 


For the term of angular momentum, Eq.(25b), we have 


(35) 


n i,k,fi 

E ? 1" fc 
e^kafi og. 

i,k,n 

Substituting these results into Eq. (26) the Hamiltonian becomes 

Hnc = ^ + z|wS6»(n)eifcab^a^ 

i,n 

or 


(36) 


(37) 


Hnc 




■9 2 


e{n) (ai 


t 2 2t 

Of, - 



(38) 


Let us now define new creation and annihilation At operators as 


-^n ~ (®n 


yllt _ 
^B — 

1 / it , ■ 2t\ 

^ (^Ob +^afl J , 

(39a) 

^ (4 + ial) 


42 t _ 
^B — 

^ (4^ - *4^) : 

(39b) 

'B = + ^b) 


it 

Ob = 

^ (^4 + 4l4) . 

(40a) 

1 = iT - Ai) 


„2t 
®B “ 

-^B ) ■ 

(40b) 


thus 
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The components of the fields and momenta, Eqs. (33), in terms of these new operators are 


= 


= 


Tri — 


= 




2^ 

i 

'i'\/ ^fi 


{aI - Al - AlJ + ^ 


i^A]^ — A\fi + A^ — A^_f^ 


( 4it 
2 V " 


2t 


A^_n-Atn-A^ 


The operators A''^ and At satisfy the following commutation relations 


Ai 4 i 1 _ [ 4 * 1 AO 1 


= 0 , 


4* 40 ' 




Substituting Eqs. (40) into Eq. (38) we obtain the Hamiltonian in terms of At' and At 


Hnc = ^ 

n 

= E 


uofitlftiAl^Al + Al^Al) - {Al^Al - Al^Ai) )■ , 


LUn 




gujfi0{n) 


Aff, Aff + ( Hfj + 


gujf;0(n)) 2f a2 


That can be written as 


ffNC = 5]c.s[A|fVi + A|lV|] 


(41a) 

(41b) 

(41c) 

(41d) 

(42) 


(43) 

(44) 


where we defined the quantum numbers A^i = Al'Al and NS 




besides the terms of energy deformation 


A^ — t/1+ 


At = 


/ guifi0{n )\ gujn0{n) 

V 2 


/i , (9^n0{n)\^ ^ guJnO{n) 

V +^— 


(45a) 

(45b) 


We can see that the free theory of the non-commutative scalar field with two real scalar components leads to the 
emergence of two types of particles, which correspond to particle and antiparticle that are no more degenerated [21]. 
There are small corrections, parametrized by 9 and cr, on their standard energy expression E = where 

k = 27r|n|/i? is the wave number. 

Taking the energy relations El = ojrAl for particles and ES = Wr:A| for antiparticles descriptions, we consider Nl 
and Ni being the number operators of bosons and antibosons with momentum k, respectively. This enables us to 
analyze the thermodynamic properties of a charged and uncharged non-commutative bosonic gas. 


III. THERMODYNAMICS OF THE NON-COMMUTATIVE FIELDS 

In the previous section we explored the implications that the non-commutativity of the target space causes in the 
dispersion relation of a system constituted by massive scalar bosons. In this section we will formulate the quantum 
statistical problem for a non-commutative gas within a volume V for the purpose of studying its thermodynamics 
investigating how the introduction of the non-commutativity changes the system. We shall focus our studies in 
obtaining relevant thermodynamic quantities such as pressure, internal energy, particle number, specific heat, and the 
Bose-Einstein condensate in the ultra-relativistic (UR) and non-relativistic (NR) limits. 
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A. Thermodynamics of massive bosons gas 


We will first consider a situation that the system has no antibosons. In grand canonical ensemble we use the 
Hamiltonian defined in Eq. (44) written in terms of the momentum k 

2 


n 


—^CJrAi /Vi I ’ 

1 — 26 ^ fc/ \1 — 26 ^ fc/ 


(46) 


where N = + -^|) i® the operator for the total number of bosons and we are using the system of units 

ks = fi = c = 1; z = is the fugacity, /r is the chemical potential and /3 = 1/T. 

Thus, the natural logarithm of the grand partition function is 


InS = - 


|^ln(l — ze + ln(l — ze 


(47) 


We are dealing with a gas formed by N harmonic oscillators within a volume V ^ B?. In the thermodynamic limit 
E —> 00 (i? —> oo) the sum is replaced by an integral in the momenta 




(Pk 


Et f dkP. 


This makes the Eq. (47) given by 


\nE{/3,V,z) = - 


V 

2'k‘^ 


dkP 


(27r)® 2P 


ln(l-ze-^^'W)+ln(l-ze-^^'«) 


where [41], 




E‘^{k) = uj\l 1 + 


je{k) 

Stt 


26»(fc) 

Stt 


with 


(48) 

(49a) 

(49b) 


UJ 


= \/fc2 + m? and 6{k) = 6e ^ 




(50) 


Given the complexity of the quantum statistical problem here formulated, it is not possible to find exact analytic 
solutions for Eq. (48), because if*(fc) is a non trivial function of k, 9 and tr. An alternative to overcome these difficulties 
is to find suitable approximations able to provide solutions that enable the advancement of our investigations with 
analytic calculations of the system properties under study. This procedure is made through ultra-relativistic and 
non-relativistic limits. Alternatively we also treat the problem exactly using numerical approaches. 

In this paper we obtain the non-commutative bosonic gas properties developing both procedures, numerical and 
approximate, in view of checking validation limits of the UR and NR approximation. The relevant quantities for this 
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investigation are the total number of bosons N, the internal energy U and the pressure p given by the equations 


N = z^JnE{/3,V,z), 
U = -^lnE{(3,V,z), 

P = 

We can also analyze the specific heat at constant volume defined by 


1 fdU 


N \dT 


N 


where the partial derivative dUjdji can be calculated through 

_ /(m\ /m\ 

) V \ ^/^ / z,V \ J /3,V 
'du\ _ /d[/\ 

. ) z,v \ J py 

Replacing the equation (48) into (51) we find the expressions 


m" 

dp/V 

(h\ 

dp) N,V 

dN\ 

dp) z,v, 


N = 


V 

27r^ 


dkk^ 


1 


dN 

dz 


P,v 


^-iQPEHk) _ I ^-iQpEHk) _ I 


(51a) 

(51b) 

(51c) 

(52) 


(53) 


(54a) 


U = 


V 

27r^ 


dkk^ 


E^{k) 


E^{k) 


^-l^PE^k) _ I 2-le/3£)2(fc) _ 2 


(54b) 


1 


P = 


27r^/3 


J dkk^ {ln(l - + ln(l - . 


(54c) 


Here, the total number of bosons N = Ni -\- N 2 is a conserved quantity and an important realization that must be 
done is to make and N'^ positive definite, so that we find the transition temperature Tq at which a relatively high 
fraction of ‘atoms’ begins to condensate on the lower energy state by taking 


, . ,'m0\ wP9 


(55) 


For temperatures greater then Tq we must have p, < m! and the system becomes completely excited. Starting from 
the equations (54) we calculate the specific heat via partial derivative dU/dp given by Eq. (53) for T > Tg by using 


dir 

dp / z,v 


V 


dkr 


E^{k) 


27r^ ./q 1 V — 1 J \ z~'^edE'^(k) _ 


^-l^PE\k) ^ 


E^{k) 




dU 

dz 




1 rdN 


^- 2 ”^ ,2 i E\k)z-rf^^"r E^{k)z-rf^^"r 

' dkr > ^ ' I 5 ^ 


z,V 


-^0 1 {z~^edE'^(k) _ ^^-lQPE'^(k) _ I'j2 P 


dN \ 

Ih ) 




— 1 roo 


p,v 


2r 


dkr 


^-1^0E\k) 


^-l^pE^ik) 


(z-le/3SMfc) - 1)2 (z-le/3i5"(fc) - 1)2 [■ 


For T < To we use 

/aC/\ V f E^(k)(E^(k)-m')z-r>^^'r E^(k) (E^(k) - m') 

[dp)y ~ 2r Jo \ - 1 )^ (^z-iQpE^k) _ lY 
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where 2 = . 


1. The ultra-relativistic limit 


The approximate calculations are made from the expression (48) by expanding the integrand as a power series of 
the fugacity z 


InS = 


V 

27r^ 


dkk 




-PE\k) 1^2 ^^-2i3Ekk) _i_^-20E^{k) 


+ ... . 


(58) 


By using this approach, the ultra-relativistic limit is computed by making m Tq 1), so we expand the 

equation above up to first order in j3m and then we expand in the parameters 9 and a. Assuming 9j< \ and 
a!j5 < 1 we find as result that the main terms to this approximation are 


lnH(/3,T,z) 


V 

27r^ 



-f 


756(2 


22056(2^2 

4/3 


53 ( 2 ) -f 


330750^ 
1024/3'^TT' 


53 ( 2 ) 


6548856I^(j2 28378350® , A 

64^9^4 9w{z)+ ^2334^9^6I • 

(59) 


Then from the partition function we find the equations below at ultra-relativistic limit 


N = 


U = 


V 

2ii^ \ /3® 

V r 12 


0353 ( 2 ) -F 


27r2 (/34 


1 


/Sj 54(2) -h 


7502 

8/3®7r' 


37502 

8/3®7r" 


7502 


2„2 


55 ( 2 ) - 


53 ( 2 ) - 


220502(7 

4/3 

1543502(7 

4/3®'7r2 


57 ( 2 ) -f 


330750"^ 

1024/3^774 


4_2 


57 ( 2 ) - 


65488504(7 

64/3®7r4 


• 59 ( 2 ) -f 


28378350® 

16384/3977® 


59 ( 2 ) 


(60a) 


(^) + 


2_2 


2772 (^ 3454 ( 2 )+ 


53 ( 2 ) - 


220502(7 

4/3®772 


53 ( 2 ) -f 


23152504 

1024/38774 

3307504 

1024 / 38774 ® 


'4_2 


53 ( 2 ) - 


589396504a 

64 / 31 O 774 

65488504a2 


510 ( 2 ) -f 


255405150® 

16384 / 32077 ® 


510 ( 2 ) 


64/31®774 


510 ( 2 ) -I- 


28378350® 

16384^2077® 


510 ( 2 ) 


(60b) 


(60c) 


where ga(z) is the Bose function defined as 

n—1 


This function diverges at 2 = 1 when a < 1, but becomes the Riemann zeta function ((a) when a > 1. 

The specific heat (52) is computed using the equations (60a) and (60b). As we previously stated for T < Tq we 
have 2 = e ^”2 then 


( 



V 


+ 


T f 12 




37502 

8/3®772 




^55405150® 
V16384/32077® 


58939650V\ 
64 / 32 O 774 J 




-f 


/ 23152504 
V 1024^8774 


1543502a2\ 
4/38779 J 




For T > To we use the following equations 
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2. Nan-relativistic limit 

To take the non-relativistic limit m 3> Tq and we can expand the modified energies (49) in powers of /vn? up to 
first order and with respect to the parameters 0 and a. We obtain the following results as the major contribution to 
NR limit 

n^e 


E\k) 

“ "‘+SS 

E^k) 

k? 



+ 


mr9 

Stt 


Now we substitute the energies above into Eq. (58), such that we find 

3/2 


In; 


E = V(^\ 


\2tt/3 J 


95/2 ze 




+ 95/2 I ze 




Given the partition function we find the expressions for the total number of particles and internal energy 

3/2 r 


N = V 


y 

\2nf3j 


93/2 


- 93/2 




, m'^9\ ( rn?9\ W [ m \ 

^ ' Sir J ^ V Stt J 2P \2ttI3J 


3/2 


95/2 2:e V / + 35/2 ( 0e V ) 


(63a) 

(63b) 

(64) 


(65) 


( 66 ) 


where we see the non-commutative terms acting differently in the rest mass of the two types of bosons. In this 
approximation, the critical chemical potential and the critical temperature are given by 


9o = m- 


'mr9 

Stt 


and N = V 




(wo. 


3/2 r 


93/2 (1) + 93/2 ( e 


,-/5o 


(67a) 


We can see that the non-commutative theory modifies the thermodynamics of the gas introducing corrections that 
does not exist in the standard BEG. All the previous non-relativistic results known in the literature are recovered as 
9^Q. 


B. Thermodynamics of massive bosons-antibosons gas 


To describe the system with antibosons we express S = Tr where Q = ~ -^|) i® ^ operator 

which corresponds to a conserved quantum number generically referred as to charge. The grand partition function is 
now explicitly given by 


00 



( 68 ) 
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Recall that we have previously identified A| and A| in (45) as the energies of the particle and antiparticle, as well 
discussed in [21]. The grand partition function in the thermodynamic limit is 




dkk‘' 


ln(l - + ln(l - 


(69) 


Analogously to the previous section, we obtain the expressions for the net charge Q, the internal energy U and the 
pressure p through the equations 


2^1nS(/3, U,2;), 

(70a) 

-^lnS(^,U,z), 

(70b) 

^ln5(/3,U,z), 

(70c) 


and the specific heat at constant volume is 


cy = 


-/32 fdU 


Q \ d/3 


Replacing the equation (69) into (70a) we find the expression for the charge 


Q = 


V 

27r^ 


dke 




(71) 


(72) 


In this context, we have a conservation of the number Q = Ni — N 2 , where A^i and N 2 are the total number of bosons 
and antibosons, respectively. Making and N'^ positive definite we get to the conclusion 


-U}\ 1 + 


jQ/k) 

Stt 


^e{k) 

Stt 


< p < UJ\ 1 + 


j9{k) 

Stt 


'‘^d(k) 
Stt ’ 


(73) 


where we suppose that the bosons number exceeds the antibosons number, i.e., Ni > N 2 and thus Q > 0, then we 
find p > —a;^0(fc)/87r and the transition temperature Tg is obtained by taking 

' 2 


k = m' =m\ll + { —j - 


m9\ m39 
Stt 


(74) 


Note that the Eq.(72) is actually the number of excited charges, Q — Qo, where Qo is the number of charge in the 
ground state. So at /3 = /3o and /r = to' we can get total charge density p = Q/V written as 

' ui‘^e(k) 
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27r^ 


sinh 
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Using Eq.(70b) we get the internal energy of the charged system 
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(75) 
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(76) 


which can also be written as 
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( 77 ) 
















































and the pressure is 
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1 


P = 


dkk^ 


ln(l - + ln(l - 


27r2/3 Jq 

Starting from the equations (72) and (76) we can calculate the specific heat. For T > Tq 

2 / / 7 \ \ 2 
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(78) 


(79a) 


(79b) 


(79c) 


and for T < Tq 
'dU 


V f°° ^,,2 i ^Hk) (E^ik) - to') E'^{k) {E'^{k) + to') ze^^'W 

dpj^ ^1 \ r^^F777^^72 + r^i;^77^^72 


p-ie/3£dfe) - 1 )^ 


[zePE^{k) _ ly 


where we have z = . 


1. The ultra-relativistic limit 


The approximate calculations are obtained from the expression (69) by expanding the integrand as a power series 
of the fugacity zi = z and Z 2 = z~^ 


InS 


V 

27r^ 



e-/5i5bfc)^^ + Ie-^13E^(>^) 



+ ie-^0E\k) 



(80) 


This is analogous to section III A 1, we take the ultra-relativistic limit and then we find 


\nE{/3,V,z) 


^{|l b4(^)+54(2 ")]+^[95{z)-g5iz [97 {z) - giiz ^)] 

[ff 6 (^) + 9,{z-^)] + [97iz) - gy(z-i)] I . 


( 81 ) 
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Again, since the partition function is obtained we can determine the excited charge density pe, internal energy density 
u = U/V and pressure p 


Pe = [gsiz) - gsiz ^)] + [giiz)+g4iz ^)] - [ge{z) + ge{z ^)] 


7502 


2/3® TT^ 
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P = [ 54 ( 2 ) + < 74(2 ')] + [< 75 ( 2 ) - 55(2 ')] - 11^ [- 77 ( 2 ) - griz ')] 


77^/34 
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(82c) 


But in the ultra-relativistic limit (3m <C 1 => (3p <C 1 where p < to', then 

e=^^''«l±/3p ^ Pa « 5a(l) ±/3/7ga-i(l). 


(83) 


This approximation leads to the thermodynamic quantities 
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So the expression for the total charge density is 
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(84a) 


(84b) 


(84c) 


(85) 


Starting from this equation we can obtain the critical temperature in terms of p. When 0 —>■ 0 we get Tg = 
which is in agreement with the Haber and Weldon result [28]. As consequence, for the to —>■ 0 limit we get that Tg —?► 00 
and hence all net charge of the massless bosons gas resides in the ground state, differently from the non-commutative 
gas where if to —>■ 0 we will have Tg ~ (7r®p/30C(4)) 
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In addition, the specific heat for T > Tq can be found by 


dp y zy 


m 

dz 

dz 


py 


py 


V 1890/i6»cr2 11256|2 11025^03 I 
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(86a) 


(86b) 


(86c) 


and for T < Tq 


dU\ V(24^,^, 48m'e^,^, 1620m'6»cr2 11256»2 4725m'6»3 , 
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2. The non-relativistic limit 


In the non-relativistic limit we use the energies (63). We subsequently substitute these energies into Eq. (80), to 
find 


InS = V 



3/2 


95/2 



+ 95/2 



(87) 


In the case of non-relativistic limit the temperature is low enough such that /3m 1, then the second term into equa¬ 

tion above corresponding to antibosons can be neglected. But still has non-commutative modification in the particle 
term acting as an additional mass which decreases the boson rest mass. In this approximation the thermodynamic 
equations are 



The specific heat reads 

_ 15 g5/2(^) _ 9 93 / 2 { z ) 
4 93/2{z) 4gi/2(z)’ 


(88a) 

(88b) 

(88c) 


(89) 


where we have set z = ze V S’" J _ z = 1 for T <To, while z < 1 for T > Tq, we fall back at the situation of usual 
non-relativistic Bose gas. 

Therefore, in the study of the Bose-Einstein condensation through non-commutative fields arise corrections in the 
system properties, which can open new investigations such as looking for small deformation effects produced by such 
systems in high energy physics, astrophysics, cosmology or even in condensed matter. 
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IV. RESULTS AND DISCUSSIONS 

The thermodynamic behaviour of scalar fields was explored in the context of quantum held theory on a non- 
commutative target space, as presented in section III. This leads to deformation of physical quantities that will 
be analyzed in detail in this section. The results are compared to those well-established in the literature by the 
usual commutative theory. We shall plot graphics of the modihed thermodynamic functions in different regions of 
temperature, which show how the non-commutative parameters B and a affect these quantities. These parameters are 
expected to be related to short distances, thus they cannot be large and were set at the following values 0 = 2.4 x 10“^^ 
and (T = 1.1 X 10“^®. For convenience we assume the value of parameters in arbitrary units. 


A. Bose-Einstein condensation withont antibosons 

We will present the contributions made on a system of iV fixed massive particles. The critical temperature Tq at 
which the Bose-Einstein condensation occurs corresponds to fj, = m'. For some temperature above Tq we can always 
find a chemical potential satisfying ^ < m' such that Eq.(54a) holds. On the other hand, for T < Tq the ground 
state becomes macroscopically occupied, then the number of particles in excited states decreases proportionally. This 
behavior can be seen in the Figure [2] for commutative and non-commutative bosons with the mass adopted as 
m = 5.1 X 10®. 



FIG. 2: Comparison of the ratio of particles in the ground state as a function of the temperature for commutative and non- 
commutative bosons. Fixing 6 = 2.4 x cr = 1.1 x 10“^® and m = 5.1 x 10® we choose n = N/V = 1.0 x lO"®’ thus 

To « 3.26 X 10^® and « 3.45 x 10®®.“ 

“We will refer Tc and Tq as the critical temperatures of commutative and non-commutative cases, respectively. 

Analyzing the Figure [2] we note that the behavior obtained after introduction of non-commutativity is similar to the 
commutative case. However, the critical temperature is modified which can also be seen in Figure [3], where are shown 
the numerical results of Tq for different values of particle density, as well as those obtained by the non-commutative 
ultra-relativistic (UR-NC) regime. For a given particle density there is a temperature region where the transition 
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temperatures of non-commutative curves are lower. This effect begins to appear in the graph in high temperatures 
(T ~ 10^^). We can also check the validity of our UR-NC approach that coincides with the non-commutative exact 
result in a limited range of temperature. 



FIG. 3: Results of To as a function of density particle n comparing commutative (exact-C), non-commutative (exact-NC) and 
ultra-relativistic (UR-NC) cases. 

Continuing our investigations we study the non-commutative effects on internal energy and specific heat, whose 
behaviors are shown in Figure [4]. The results show an increase in the discontinuity dX T — Tq of specific heat of 
the system, that features a second-order phase transition. The gap is equal to 8.55 for non-commutative bosons 
while is 6.26 for commutative bosons. Moreover, it is interesting to observe that in the region T > Tq arise a lump. 
A contribution that can be related to the angular momentum term in the non-commutative Hamiltonian. In low 
temperature physics, these effects can be even more interesting in condensed matter physics, such as in high-Tc 
superconductivity and thermodynamic anomaly of water. 


B. Bose-Einstein condensation with antibosons 

Taking into account the possibility of boson-antiboson pair creation, as presented in the section IIIB, we fixed the 
charge Q so that the transition temperature Tq can be obtained numerically by equation (75). The results for ratio 
of charges in the ground state as a function of the temperature is presented in Figure [5], where one can be seen that 
non-commutative effects appears at lower temperatures than for the system without antibosons. It is interesting the 
way the critical temperature varies according to the choice of the parameters 6 and a. An increase in 6 causes a 
decrease in the transition temperature, but when 9 approaches zero or a is large enough the non-commutative curve 
coincides with the commutative curve what is in accord to the regularization given by a Gaussian-type distribution 
in Eq. (17). 

The results of Tq for different values of charge density can be seen in Figure [6] comparing the numerical procedure 
with UR-NC regime, which shows the efficiency region of this approximation based on the expansion of the parameters 
9 and a in determining the thermodynamic quantities of interest. Another interesting aspect presented by this new 
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T/To 

(a) 



(b) 

FIG. 4: (a) Internal energy per particle as a function of the temperature, (b) Specific heat as a function of temperature with 
a fixed number of particles. Using the same parameters of Fig. [2]. 

non-commutative BEC arises in the curves of transition temperatures versus charge density for different values of 
mass shown in Figure [7]. 

The non-commutative effects on internal energy, pressure and specific heat are shown in Figure [8] together with 
those obtained through the usual commutative theory. In this regime of temperature, the numerical results obtained 
using the equations (76), (78) and (79) are the same as those obtained by the ultra-relativistic regime. We can observe 
that the specific heat has no gap for a fixed net charge Q. This is because the mechanism of pair creation of particles 
and antiparticles removes this anomaly in three spatial dimensions [36]. 
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(a) 




(b) (c) 

FIG. 5: Results for the ratio of charges in the ground state as a function of the temperature choosing p = 6.12 x 10^"^ and 
m = 5.1 X 10®. (a) Comparing commutative and non-commutative cases, the critical temperatures are Tc ~ 6.0 x 10® for 
commutative bosons and To ~ 5.1 x 10® for non-commutative bosons with 0 = 2.4 x and cr = 1.1 x 10~^®. (b) Setting 

a = 1.1 X 10“^® and varying 6. The dotted curve presents a critical temperature To « 5.85 x 10® while the dashed line 
presents Tq ~ 3.55 x 10®. (b) Setting 6 = 2.4 x 10“^^ and varying a. The solid, dashed and dotted lines have Tq « 5.1 x 10®, 
To ~ 5.46 X 10® and To ~ 5.84 x 10®, respectively. 
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FIG. 6: Numeric results of To as a function of charge density p, comparing commutative (exact-C), non-commutative (exact-NC) 
and non-commutative ultra-relativistic (UR-NC) cases. The parameters are 9 — 2.4 x a = 1.1 x 10“^® and m = 5.1 x 10®. 



FIG. 7: Numeric results of To as a function of charge density p for bosons with various values of mass m = 5.1 x 10^, 
m = 5.1 X 10®, m = 5.1 x 10®. The solids lines represent the non-commutative situation and the dashed lines indicate the 


commutative case. 









(sOx) °xO/n 


24 



(a) (b) 



(c) 


FIG. 8: Comparing the commutative and non-commutative cases with a net fixed charge density p = 6.12 x 10^^, setting 
6 — 2.4 X 10“^^ and a = 1.1 x 10“^®. (a) Internal energy per particle as a function of the temperature, (b) Pressure per density 
of particles as a function of the temperature, (c) Specific heat as a function of temperature. 
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V. CONCLUSIONS 

We investigated the thermodynamic properties of a relativistic boson gas in the scenario of quantum field theory 
based on non-commutative fields. We aimed to obtain the thermodynamic equations that describe a Bose-Einstein 
condensate, such as temperature transition, fraction condensate, internal energy, pressure and specific heat, in order 
to investigate the corrections caused in the system due to the introduction of the non-commutativity. 

The first consequence of the non-commutativity appears in the dispersion relation, which leads to modifications 
of the thermodynamic quantities. To understand these modifications and their dependence on the non-commutative 
parameters {9 and cr) introduced in the theory we mainly study the graphics obtained numerically and by approximated 
methods at some regimes. These modified thermodynamic quantities are then compared to known results of the 
commutative theory in the literature. Specially, in the case without antibosons, the specific heat around the critical 
temperature has a larger gap in the non-commutative case. In low temperature physics this effect can be even more 
interesting in high-Tc superconductors. The specific heat also presents a lump for temperatures immediately larger 
than the critical temperature which are very similar to the thermodynamic anomaly of water. 

The analysis of these results shows that the non-commutative effects cause significant contributions in high critical 
temperatures, what is related to the fact that the non-commutativity is an effect only perceptible at the Planck scale, 
because the non-commutative parameters are expected to be very small in the context of quantized spacetime. The 
system with antibosons has non-commutative effects at relatively lower critical temperatures compared to the system 
without antibosons. In this case the specific heat has no gap for a fixed charge Q. This is because the mechanism 
of pair creation of particles removes this anomaly in three spatial dimensions. In this case one can recover the well- 
known Haber and Weldon critical temperature in the limit as 0 —)■ 0. On the other hand, in the limit of massless 
non-commutative scalar fields one can still find a finite critical temperature depending on 6 for the non-commutative 
Bose-Einstein condensate. 

Among the prospects to apply this work is to explore the Bose-Einstein condensation with non-commutative scalar 
fields in curved spaces. 
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Appendix A: The commutation relations 

In this appendix we calculate the commutation relations that appear in the section (IIB 2). Using the equations 
(15) we found 






- ie“6'7rc(f,t),(/?'’(y,t) - ^e“6»7rd(y, t) 

[(^“(f, t), t)] - ie“6»[(/j“(x,t),7rd(j/,t)] 
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it’^96{x-y), 

- ie“6'7rc(f,t),7r6(y,t) 

[(/?“(£, t), TTbiy, t)] - [7rc(x, t), nb{y, t)] 
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[TTaix,t),TTbiy,t)] = [^alx, t) , TTbiy ,t)] = 0. 

The commutation relations components for the Fourier transformations are obtained using (20) 
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To obtain the commutation relations for the operators ah' and ah we use equations (30a) and (30b) 
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In the following we calculate the commutation relations (42) using the expressions (39), then 
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Therefore, 
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